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ABSTRACT 



Aims. We determine the steady-state of an axisymmetric thin accretion disc with an internal dynamo around a magnetised star. 
Methods. Starting from the vertically integrated equations of magnetohydrodynamics we derive a single ordinary differential equation 
for a thin accretion disc around a massive magnetic dipole and integrate this equation numerically from the outside inwards. 
Results. Our numerical solution shows that the torque between the star and the accretion disc is dominated by the contribution from 
the dynamo in the disc. The location of the inner edge of the accretion disc varies between R A and 10/? a depending mainly on the 
strength and direction of the magnetic field generated by the dynamo in the disc 
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1. Introduction 

In this paper we present a new solution for an accretion disc 
around a magnetic star. This star could be a neutron star, a white 
dwarf, or a T Tauri-star, but we assume that it is a neutron star 
since it is easy to measure the torque between the neutron star 
and the accretion disc by timing the X-ray pulses from the neu- 
tron star. The new feature of our solution is to include the effect 
of an internal dynamo in the accretion disc. By doing this we 
hope to be able to explain the torque reversals that have been 
observed in some X-ray pulsars. 



Shakura & Sunyaev (1973)| formulated the standard model 
of a geometrically thin, optically thick accretion disc. They 
were able to obtain an analytical solution of the height- 
integrated hydrodynamic equations, after having introduced the 
a-prescription for the turbulent stress, which transports the an- 
gular momentum outwards through the disc; however, they did 
not explain why the disc is turbulent in the first place, since a 
disc in Keplerian rotation is stable according to Rayleigh's cri- 
terion. Balbus and Hawley (1991) instead showed that it is un- 



stable if there is a weak magnetic field in the disc. Subsequent 



numerical simulations (e.g. Hawley, Gammie & Balbus 1995 
Balb us and Hawley 1998[ ) confirmed that this instability gener- 



ates turbulence and that the resulting turbulent stresses transport 
angular momentum outwards. 

The interaction between a magnetised star and a surround- 
ing accretion disc is one of the most poorly understood aspects 
of accretion. At the same time, it is central for our understand- 
ing of the spin evolution of objects as diverse as T Tauri stars and 
X-ray pulsars. The magnetic field of the star penetrates the sur- 
rounding accretion disc and couples the two. According to the 
Ghosh & Lamb ( 119791 ) model, the part of the accretion disc that 
is located inside the corotation radius provides a spin up torque 
on the star, since it is rotating faster than the star, while the more 
slowly rotating outer part of the accretion disc brakes the star. 
The net torque is determined by the location of the inner edge 



of the disc, which moves inwards as the accretion rate increases, 
thereby increasing the spin up-torque on the star. 

Campbell (1997) proposes physical descriptions of the mag- 
netic diffusivity in terms of turbulence or buoyancy, and 



Campbell & Heptinstall (1998) solve the resulting equations nu- 
merically. For both forms of diffusivity, the magnetic coupling 
between the disc and the star leads to an enhanced dissipation 
in the inner part of the accretion disc compared to the stan- 
dard |Shakura & Sunyaev (1973)| model. This raises the tempera- 
ture such that electron scattering dominates Kramer's opacity at 
larger radii than is otherwise the case, thus increasing the frac- 
tion of the disc that is subject to the (Lightman & Eardley 1974) 
instability. Brandenburg & Campbell (1998)| considered a form 
of magnetic diffusivity that allows further analytical progress to 
be made, but the qualitative results remain the same. 

All the models predict a positive correlation between the ac- 
cretion rate and the torque on the neutron star and even pre- 
dict a negative torque on the neutron star at very low accretion 
rates. Timing of X-ray pulsars during outbursts of Be/X-ray tran- 
sients have provided at least some qualitative support for such 
a correlation (e.g. Par mar et al. 19 89). The BATSE instrument 
on the Compton Gamma Ray Observatory made it possible to 
extend this database significantly (Bildsten et al. 19971 . In par- 
ticular there are a few X-ray pulsars with permanent discs that 
are oscillating between phases of constant spin-up and constant 
spin-down without a significant difference in the X-ray luminos- 
ity between these states, which appears to contradict the standard 
model for a disc-accreting X-ray pulsar. 

|Nelson et al. (1997)] propose that these torque -reversals can 
be the result of transitions between co-rotating and counter- 
rotating accretion discs, though several other models have also 
been proposed. Torkelsson (1998) argue that the torque between 
an accreting star and its disc can be enhanced by the presence 
of a magnetic field generated by the turbulence in the accretion 
disc. The torque reversals are then the result of a reversal of the 
magnetic field generated by this dynamo. However, he did not 
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construct a self-consistent model of the accretion disc. The aim 
of this paper is to construct such a model of an accretion disc 
with an internal dynamo around a magnetic star. We work in the 
spirit of Shakura & Sun yaev (1973)| and assume that the disc is 
geometrically thin. In Sect. 2 we start from the equations of mag- 
netohydrodynamics (MHD) and derive a single ordinary differ- 
ential equation for the radial structure of the accretion disc. We 
then present numerical solutions of this equation in Sect. 3 and 
discuss the properties of these solutions in Sect. 4. Finally we 
summarise our conclusions in Sect. 5. 



2. Mathematical formulation 

We study a steady, thin axisymmetric Keplerian disc around a 
star with a magnetic dipole field. The basic equations describing 
the structure of the thin accretion disc can be derived from the 
equations of magnetohydrodynamics. 

2. 1 . Conservation of mass 

In steady state the continuity equation takes the form 

V-(pv) = (1) 

where p is the density and v = (Vr, v^, v z ) is the fluid velocity 
with radial, azimuthal, and vertical components , respectively. 
For a thin axisymmetric disc and after neglecting a vertical out- 
flow from the disc, we get 



1 <9 

-— (RT.v R ) = 0, 
RdR 

where 2 is the surface density 
2= f pdz^2pH, 

J-H 



(2) 



(3) 



and H is the halfthickness of the disc. For a steady disc the inte- 
gral of Eq.© gives 



r" 

M = -2nR I pv R dz = -2nRl,v R = constant, 

J-H 

which is the accretion rate. 

2.2. Conservation of momentum 



(4) 



If assuming a steady state the Navier-Stoke's equation can be 
written as 



p(v V)v = -VP+pVO + JxB 



+V-|pv|Vv- -(V-v) 



(5) 



where P is pressure, v kinematic viscosity, cD the gravitational 
potential J = — (V x B) = (J R , J^, J z ) the current density, and 
B = (B R , Bq, B z ) the magnetic field. The viscosity is in general 
low, and we only retain it where it plays a crucial role. 
The radial component of Navier- Stoke's equation is 



Vr 



dv R v <f, 
dR R 



Ho \ dz 



8B, 
dR 



dP pGMR 
~dR ~ (R 2 +z 2 ) 3 ' 2 ' 



(6) 



For a thin accretion disc » c s as shown below and the domi- 
nant terms of the equation give us 



*+ R ' 



(7) 



which shows that the disc rotates in a Keplerian fashion. 

In similar manner, the vertical component of the momentum 
equation for a steady flow is 



dv z dv, 
dR z dz 



B,/, dB<f, 
Ho dz 



B R dB R B R dB z 
Ho dz ho dR 



dP 
' dz 



pGM 



(8) 



Neglecting vertical outflows and assuming the magnetic field to 
be weak the equation reduces to the equation of hydrostatic equi- 
librium 



1 dP 

p dz 



GM z 
' R 2 R' 



(9) 



Using H as the halfthickness of the disc, the pressure at the mid- 
plane of the disc is 



2 R 3 

but the hydrostatic equilibrium can also be expressed as 



H 
~R 



Cs 
Vkepl 



(10) 



(ID 



which shows that the Keplerian velocity is highly supersonic in 
a thin accretion disc, as assumed above. 

The azimuthal component of Navier-Stoke's equation re- 
duces to 



vj_d_ 
R dR 



Br18_ 
/to RdR* 



B z 8B$ 
Ho dz 



+ R 2 dR 



RpV dRW 



(12) 



We neglect ^■■^(RB ( p) because the the radial length scale is 
much longer than the vertical length scale in a thin accretion 
disc. Integrating Eq. (fT2l) vertically across the disc and multiply- 
ing both sides by R, we get 



B 7 Ba 



A'o 



1 d 
R+- — 

RdR 



dR \R 2 



(13) 



where the specific angular momentum I = Rv$ oc R x l 2 . The 
magnetic term describes the exchange of angular momentum be- 
tween the disc and the star via the magnetosphere. This term 
vanishes if B^ is an even function of z, but the shear between the 
disc and the stellar magnetosphere generates an odd B^ whose 
value in the upper half of the disc is 



B,i, 



= -yB 



(14) 



where Q.^ = v^/R, O s is the angular velocity of the star, and y is 
a dimensionless parameter of a few (Ghosh & Lamb [l979t . 

In this paper we consider the effect of adding a large-scale 
toroidal field that is generated by an internal dynamo in the 
accretion disc. Such a dynamo is a natural consequence of 
the magnetohydrodynamic turbulence in the accretion disc (e.g. 
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Balb us and Hawley 1998[ ). To estimate the size of B^a, we as- 
sume for the moment that the viscous stress in the accretion disc 
is due to the internal magnetic stress 



fM) 



= a ss P, 



(15) 



Shakura & Sunyaev (1973)| prescription for the 



where we use the 

viscosity in the last equality. Based on the results of numeri- 
cal simulations of magnetohydrodynamic turbulence in accre- 
tion discs (e.g. Brandenburg et al. 1995) Torkelsson (1998) ar- 
gues that 



Tdyn = 



B4, 



B,p 
~B~ 



(16) 



where ydyn ~ 10. However, this B^ is the sum of the large-scale 
field and a small-scale turbulent field, that is also contributing 
to the stress through its correlation with a turbulent Bs-field. 
Since the large-scale field might be a small fraction of the total 
field we multiply B$ with a factor e to get an estimate for B^dyn : 



#0,dyn = £ (assAWdyn^) 



1/2 



(17) 



where -1 < e < 1, and a negative value describes a magnetic 
field which is pointing in the negative ^-direction at the upper 
disc surface. 

We can now estimate the magnetic pressure in the accretion 
disc, to which the toroidal magnetic field is the main contributor. 
According to Eq. ( fTTl l. the pressure is approximately 



0,dyn 

2yu 



1 



(18) 



Since we use |e| < 1, a ss = 0.01 and ydyn = 10 in our models, 
we see that the magnetic pressure will not significantly affect the 
vertical structure of the accretion disc. 

The vertical magnetic field can likewise be split up into two 
components; (i) the stellar dipolar magnetic field, whose value 
in the stellar equatorial plane is 



^z.dipole - 3 , 



~R 3 



(19) 



where fi is the magnetic dipole moment; and (ii) a dynamo com- 
ponent Z^dyn- There is no obvious way to model B^dyn within our 
one-dimensional model, but numerical simulations like those by 
Brandenburg et al. ( 1995) suggest that B z and Br are compara- 
ble, SO We expect that B^dyn ~ fi^dyn/Tdyn- l n f aCt > ^<",dyn Will 

modify the structure of the poloidal magnetic field, but this can 
happen even if there is no internal dynamo in the disc because of 
the currents that are induced in the disc (e.g. Bardou & Heyvaerts 
[19961 

We can now expand the product B-B^ as 

B Z B$ - (fi z ,dipole + Bz.dyn) (^0,shear + #0,dyn) = 



^z,dipole^0,shear ^z,dipole^0,dyn 

B z ,dynB ( p f shear ^s.dyn-^^.dyn- 



(20) 



As we see below the term -B z ,dipoie^,dyn is significantly greater 
than B;,dip O iefi0,shear almost everywhere in the disc. The term 
5 z ,d yn fi0, S hear is smaller in size than B z4ipo \ e B^ dyn by a factor 
y/jdyn, and since it is difficult to model it, we ignore it. Finally 
the term B, ; d yn fi0,dyn does not contribute directly to the exchange 
of angular momentum between the accretion disc and the accre- 
tor, though it does affect it indirectly by contributing to the ra- 
dial transport of angular momentum through the disc. Although 



this term can be important, we have decided to ignore it since 
there is no obvious way to model it in our one-dimensional ap- 
proach. One should notice here that this does not change the 
main qualitative conclusion of this paper that an internal dynamo 
in the accretion disc makes a significant contribution to the ex- 
change of angular momentum between the disc and the star; on 
the contrary, we would see a stronger effect if we were to keep 
the B 7 dyn#0 dyn-term. We now expand the ZLZ^/ji/o term in Eq. 
03 



1 d 
RdR 



B z O k - Q s 
d / / 



dR \R 2 



(21) 



2.3. Conservation of energy 



For a slow inflow of matter through an optically thick disc, the 
local viscous dissipation v • f v is balanced by the radiative losses 
V • F ra d. This gives us 



9 GM \o-T A c 
8 fl 3 3r 



(22) 



where T c is the temperature at the midplane of the disc, and <x 
the Stefan-Boltzmann constant. The optical depth of the disc is 
given by 



1 

T = pHK R = ^K R , 



(23) 



and we assume that the opacity is given by Kramer's law 

K R = K 0P T; 1/2 m 2 kg- 1 , (24) 

where k = 5 x 10 2() m 5 kg- 2 Kr 7/2 . 

Equation d22l should also contain a term describing the mag- 
netic dissipation 



J 2 2B 2 
o- crp 2 } 



2e 2 a ss y Ayn riP 
H ' 



(25) 



where rj - 1/crfiQ is the magnetic diffusivity, and we have used 
Eq. (fT7b in the last equality. This should be compared to 



9 v£ GM _ 9 vP 
8 V ~W ~ 4~H' 



(26) 



so we see that the magnetic dissipation is negligible as long as 
77 * v. 

2.4. Structure equations 

We reduce these equations to a single ordinary differential equa- 
tion for the radial structure of the accretion disc. First we assume 
the equation of state of an ideal gas, 



P = 



m p fi 



(27) 



where £b is the Boltzmann constant, fi the mean molecular 
weight, and m p the mass of a proton, but the pressure can also be 
expressed using the equation of hydrostatic equilibrium 



1 I.GMH 

2 R 3 



(28) 
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which gives us a relation between H and T c 



H = 



k B \ 



1/2 



tpGM 



l/2r>3/2 



T l J L R 



The viscous stress tensor gives us the equation 



H~ l = a SSJ P(r), 



4~' \ R? ) 
which we solve for the density of the gas 



3 -,/'« P M 3/2 
P = 4^ U"J 



T 3/2 R3 



GM. 



The optical depth of the disc is 
t = K Q —a„ - — 



(GM) 3/2 (vS) / r;° J R 



16 ss \k B ) 



-6r>-9/2 



Using Eqs. ( l22l and ((32j we get 



1 / 4 Qr _1/5 ii' 1/4 - 3/10 »- 3/4 



where y = vZ and 



/243*b\ 1/10 /Gm 



C = 



\512cr 



1/4 



The pressure is then given by 

P = C x i? l% a- 9 J w M ll% y ll l™R- ni \ 
where 



3 cl/2 /243. r 1/2 °/ 



d = -G 



\512cr 



Giw, 



3/8 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



The ^-component of the magnetic field generated by the internal 
dynamo can be expressed using Eqs. (TTTb and (l35l ) as 



fi^yn = c 2e r^«L /2 V 3/I6 M 7/16 y 7/40 /?- 2I/16 , 

where 

C 2 = Gu Ci) 1/2 . 

The magnetic field due to the shear can be written as: 
fry 



^ip, shear _g 



3/2 



where 



GMP 2 . ^ 

spin 

4tt 2 



1/3 



=; 1.5 X 10 6 P 2/3 M, 1/3 m 

spin 1 



(37) 



(38) 



(39) 



(40) 



is the corotation radius, at which the Keplerian angular velocity 
is the same as the stellar angular velocity. Here f sp i n = 27r/f2 s is 
the spin period of the star, and M\ - M/M Q . Equations ( 137) 1 and 
( 1391 show that the magnetic field generated by the internal dy- 
namo varies more slowly with radius ~ R l 3 than the magnetic 
field due to shear, ~ R 3 . Thus the dynamo component domi- 
nates at large radii. 

Equation (f2Tb gives us an ordinary differential equation for y 



y 



L _ C*v lV40 R- 45nb 

6ttR 2R y 



-C A R 



-9/2 



1-1* 

Rc 



3/2 



(41) 



Table 1. Spin parameters of the models 



P(s) 


R c (m) 


w s 


7 


6.1 x 10 6 


0.71 


18.7 


1.2 x 10 7 


0.26 


100 


3.6 x 10 7 


0.05 



where 
and 



3/ioVG 
4y 



(42) 



(43) 



3/ioVG 

The solution of Eq. (|4"TT > approaches the Shakura-Sunyaev solu 
tion at large radii, thus giving us the boundary condition y — > j 
as R — > oo. 

We introduce the dimensionless variable A through 



y = AM 

and a dimensionless radial coordinate through 
R = rR A , 

where R\ is the Alfven radius, 

,2,,4 y/ 7 



Ra = 



2ttV 



GMM 2 /! 2 ) 



= 5.1 X 10 6 M7 2/7 MT 1,7 /4o 7 rn' 



(44) 



(45) 



(46) 



which is given by putting the magnetic pressure equal to the ram 
pressure of the accreting fluid (e.g. Frank, King & Raine 2002). 
Here Mb represents the mass transfer rate in units of 10 1J kg s , 
and /i2o is the stellar magnetic dipole moment in units of 
10 20 T m 3 . The boundary condition is then A — > ^ as r — > oo, 
and Eq. (l4TT > can be written as 



A 
Yr 



1 

6nr 



A' = ~ + ^- - CsA 11 ' 40 ^ 5 ' 16 - C 6 r- 9 ' 2 (1 - u s r 3 ' 2 ), (47) 



where 

c 6 = o.3 r 

and 

Ra 



4.9ea!/ 20 r f n /i 3/I6 M 1 1I/56 M- 2/35 "- 1/28 



^20 



(48) 



3/2 



ox - | y | - 6.3 M 



-3/7^-5/7 „6/7 



(49) 



is the fastness parameter. 

3. Numerical solution 

3.1. Global solutions 

We integrate Eq. (|47| | inwards from a large radius, usu- 
ally 100/?A> to which we impose the boundary condi- 
tion that A = ^, which is similar to the approach by 



Brandenburg & Campbell (1998) but most studies of accretion 



discs have rather applied a boundary condition at the inner edge 
of the disc, and then integrated the equations outwards. 

There are two possible boundary conditions that can be ap- 
plied at the inner edge of the accretion disc, either that A = 0, 
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which corresponds to that p and T 
D), or that 



d / 3 dQ 
dr \ dr 



at the inner edge (case 



(50) 



which means that the viscosity does not at all contribute to driv- 
ing the accretion at this radius. Case D is the boundary con- 
dition that has been most widely used and that was adopted 
for instance by Shakura & Sunyaev (1973)1 In this case the den- 



sity drops to zero because the inflow velocity becomes infinite, 
which is of course not realistic, and a more accurate treatment 
shows that the inflow becomes transonic close to this position 



(Paczyhski & Bisnovatyi-Kogan 1981) In case V the inflow at 



the inner edge of the accretion disc is driven completely by the 
transfer of excess angular momentum from the accreting matter 
to the stellar magnetic field (e.g. Wang 1995 1. 

For our fiducial model, we take a neutron star of M — 1 AM Q 
and a magnetic moment of 10 20 Tm 3 , which is accreting at 
10 kgs -1 . The dimensionless parameters y and ydyn are, re- 
spectively, 1 and 10 in our fiducial model, while a ss = 0.01. 
The exact values of a ss and yd yn are unimportant, since we vary 
the parameter e below, but y influences the solutions in its own 
way as is shown at the end of this section. We consider three dif- 
ferent spin periods with corresponding corotation radii and fast- 
ness parameters (see Table. 1). The system goes into the propeller 
regime for cj s > 1 (e.g. Illarionov & Sunyaev 1975 Ghosh & 
Lamb 1979). 

Firstly, for a spin period of 7 s and e = 1, 0.1, 0, —0.1, and 
-1, we get the solutions shown from the top to the bottom of Fig. 
[TJ The e = — 1 and -0.1 solutions have case D inner boundaries 
at 4.7 and 1 .6 Ra, respectively, while the other three solutions 
have case V inner boundaries (Fig. |5J. The inner boundary is 
close to R\ if e = 0, which corresponds to the absence of an in- 
ternal disc dynamo, but moves outwards as |e] increases. In case 
V, the solution continues inside the inner edge of the accretion 
disc, but the viscosity counteracts the accretion, which is instead 
driven by the magnetic stresses. This regime has been described 
by |Campbell (1998)| who discusses how the disc is disrupted 
in this region, we assume that this region belongs to a bound- 
ary layer that we do not attempt to model in this paper. On the 
other hand, all solutions approach the Shakura-Sunyaev solution 
at large radii, as required by our boundary condition. 



0.07 rr 
0.06 : 
0.05 ; 
< 0.04 

H3 °- 03 : 

0.02 ; 

0.01 : 
0.00 <- L 







1(1 



15 



20 



25 



30 



Fig. 2. jjj ( ■ft A), as a function of r for the fiducial neutron star 
with a spin peroid of 7 s and with e = 1, 0.1, and from the 
bottom to the top. 

By increasing the spin peroid to 100 s, we see that the e = 
-0.1 solution also develops a case V inner boundary (Fig. [3}, 
but there are only small quantitative changes for the e = — 1 
and -1 solutions. (The inner boundary of the e = -1 solution is 
now located at 4.4RaO The dependence of the solution on the 
spin period can be better seen in Fig. [4] where we vary the spin 
period, when e is fixed to 0. Then A has a local minimum for 
P s — 7 s, but this minimum weakens and disappears as the spin 
period is increased. This is similar to how A depends on e. For 
e = — 1, A has an unphysical negative local minima at a small 
radius, such that the physical solution has a case D boundary 
at several Alfven radii. As e is increased, the local minimum 
grows and the solution develops a case V inner boundary when 
the minimum becomes positive. Increasing e removes the local 
minimum further and the solution becomes a strictly decreasing 
function of r, which asymptotically approaches 1 /37r as required 
by our boundary condition. Increasing y has a similar effect to 
decreasing e (Fig0. 



0.20 



0.15 



0.20 prr 



0.15 - 



< 0.10 - I 




0.05 - 




o.oo 

5 10 15 20 25 30 

r 

Fig. 1. A(r) for our fiducial neutron star with a spin period of 7 s 
and e = 1,0.1,0, -0.1,-1 from the top to the bottom. 



< o.io - 



0.05 - 




o.oo -' 

5 10 15 20 25 30 

r 

Fig. 3. A(r) for our fiducial neutron star with a spin period of 
100 s and e = 1, 0.1, 0, -0.1, -1 from the top to the bottom. 

3.2. The structure of an accretion disc 

The physical structure of the accretion disc can now be expressed 
using A and r = R/Ra- We have that 

S = 3.7 x 10V /5 Mf I4 M 3 3 2/3 V 2 o /7 A(r)^ I0 r- 3 / 4 kg iiT 2 (51) 
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1 2 3 4 5 



Fig. 4. A(r) for our fiducial neutron star with e — and spin 
periods of 100, 18.7, and 7 s from the top to the bottom. 
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Fig. 5. A(r) for our fiducial neutron star with spin period of 7 s 
and e = 0.1 and y = 1, 3, and 5 ( upper dashed, dotted, and solid 
lines, respectively), e - -0.1 and y = 1, 3, and 5 ( lower dashed, 
dotted, and solid lines, respectively). 
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4000 - 
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Fig. 6. S as a function of r for the fiducial neutron star with a 
spin period of 7 s and e = - 1 (solid line), e = (dashed line), 
and e = 1 (dotted line). 
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' 1 ' " 
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Fig. 7. Vr as a function of r for the fiducial neutron star with a 
spin period of 7 s and e = — 1 (solid line), e = (dashed line), 
and e = 1 (dotted line). 



T c = 4.8 x l^V^^^^ACr) 3 ^ ^ 4 K (52) 
| = 0.012a-/ /1 V 3/8 M7 11/28 M 4 3 /35 < 4A C r ) 3/20 '- 1/8 (53) 

P = 3.o x io-V7 / V^^^ 8/3 W /14A w 11/a0r " 15/ " 

kgnT 3 (54) 

t = 370a ss 4/5 y uM 1 1 f A(r) 1/5 (55) 

v = 2.7 x 10Vf /i- 3/ X 5/14 ^f< A ('-) 3/10r3/4 m ' s_1 (56) 
v R = 84^/ 5 //- 3 / 4 M- 3/14 ^i I 3 /3 V 2( ! /7 A('-)- 7/10 ^ 1/4 m s" 1 (57) 
B,,dyn = 126 r ^a s 1 / 2 V 3/16 Mf 8 M 4 f^ 20 3/4 A(r) 17 / 40 r- 21 / 16 T (58) 

B^shear = 0.75yM 3/7 ^V 2( f V 3 (l - w. s r 3 / 2 ) T. (59) 

In Fig. [6] we show the surface density as a function of radius 
for the fiducial neutron star with a spin peroid of 7 s and e = — 1, 
0, and 1, respectively. For e = -1 the surface density attains a 
local maximum, while the other models have surface densities 



that are strictly decreasing functions of r. Figure [7] shows the 
radial velocity as a function of r. Since vr oc it becomes 
infinite as the inner edge of the disc for e = — 1, but it stays finite 
for e = and 1 and eventually goes to in the boundary layer. 
Increasing the spin period of the neutron star has a very marginal 
effect on the accretion disc, though T c , which is proportional to 
Z 2 , increases somewhat (FigJH]). 

We plot the magnetic fields for our fiducial model with vary- 
ing e and spin periods of 7 and 100 s, respectively, in Figs.|9]and 
[TOl respectively. The corotation radius, at which B,/, s hear changes 
sign, occurs inside of the inner edge of the accretion disc for 
e = 1 and -1, and for these values of e, B^dyn is the dominant 
magnetic field everywhere inside the disc. 

In general the ratio of the magnetic field due to shear to that 
generated by the internal dynamo is 

= 6.25 x 10-V^/ 2 a -i/ V 3/1 X 11/56 ^if 

^(l-ai^Air)- 17 ^- 27 ' 16 . (60) 

It is only during rather extreme conditions that the shear-induced 
field can dominate at small radii, and the dynamo is always dom- 
inant at large radii. 
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Fig. 8. T c as a function of r. The two lower curves show discs 
with e = - 1 around neutron stars with spin periods of 7 s (dashed 
line) and 100 s (solid line), respectively. The two upper curves 
show discs with e — 1 around neutron stars with spin periods of 
7 s (dashed line) and 100 s (dotted line), respectively. 
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Fig. 9. The toroidal magnetic field for the fiducial neutron star 
with a spin period of 7 s. The solid and dashed lines show the 
field generated by the dynamo for e = — 1 and 1, respectively, 
while the dotted line shows the magnetic field generated by the 
shear 
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Fig. 10. The toroidal magnetic field for the fiducial neutron star 
with a spin period of 100 s. The solid and dashed lines show the 
field generated by the dynamo for e = — 1 and 1, respectively, 
while the dotted line shows the magnetic field generated by the 
shear 



4. Discussion 

4.1. The inner edge of the accretion disc 

We have shown that there are two different forms of inner disc 
boundaries that can be found among our solutions. We denote 
these as case D and V, respectively. Case D occurs only for suf- 
ficiently negative e, while case V covers a larger part of the pa- 
rameter domain that we have studied. We summarise our values 
for the inner radius of the accretion disc in Tab. [2] We see that 
the inner radius can be significantly larger than the Alfven ra- 
dius when |e| is close to unity. This is the result of the dynamo 
enhancing the coupling of the stellar magnetic field to the accre- 
tion flow, such that it dominates the viscous torque at larger radii 
than would otherwise be the case. 

4.2. The angular momentum balance 

To understand the exchange of angular momentum between the 
accretion disc and its environment we multiply Eq. ( |2"TT i by 2ttR 
and integrate it from Rq, the inner radius of the disc, to Ri the 
outer edge of the disc 



f 



-M[^GMR x - ^GMRq) = 



4tt- 



BB, 



*' dyn R 2 dR + 



r.Rj 

Jr 



Att- 



B.B, 



—R 2 dR 



So A*o J« ^0 

-3tt (v2) Si y/GMRi + 3n (vE)« y/GMR Q . 



(61) 



The lefthand side is the difference between the angular mo- 
mentum that is advected out of the inner edge of the accretion 
disc and that which is fed into the disc at its outer edge and 
the righthand side describes the contribution of magnetic and 
viscous torques to the angular momentum balance. The term 
37t(v2)r 1 y/GMRi describes the viscous tension at the outer edge 
of the disc, which is not considered further in this paper. Rather, 
we consider the exchange of angular momentum between the 
accretion disc and the neutron star. Firstly, we have the angular 
momentum which is advected from the disc to the star 



N a dv = M y/GMRo 



2.6 x 10 2 V 20 



™l lY1 l3 ' 



(62) 



Secondly, we have the magnetic torques on the neutron star, 
which we divide into one part due to the shear 



shr 



r 

Jr, 



R < B-B, 
An— 

Ro A> 



— R 2 dR 



-7.5 x 10 26 wM /7 M?/ 7 



•J i'u 

and 

--/ 

Jri 



W20 



J/2 



)]dr 



(63) 



and a second part due to the dynamo 

r Ri B 7 B rb 

Ndyn 



An- 



^0,dyn 



R 2 dR = 



1.2 x 10 



28„ ,1/2, 



f 



A 17/40 r -37/16 rfr _ 



(64) 



Finally the viscous stress at the inner edge of the accretion disc 
transports angular momentum outwards away from the neutron 
star resulting in a torque 



-3tt (vS) Ro 

_, 7 . 
^20 



y/GMRo 



-l.A^l^^MfM^K?! 2 . 



(65) 
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This torque vanishes for a case D inner boundary. We can now 
calculate the torques on our fiducial neutron star for our choices 
of spin periods and e. These results are summarised in Table [2] 
Unless e — we see that the magnetic torque due to the dynamo 
is always significantly stronger than the magnetic torque due to 
the shear, and both are stronger for e = 0.1 than for e = 1. The 
reason for this effect is that the central hole in the disc grows 
too large when e = 1 . The dominant torque at e — lis therefore 
the viscous torque at the inner boundary, which has usually been 
ignored. 

Let us now compare our results with the BATSE data 
dBildsten et al. 19971 1 . The 7.6 s X-ray pulsar 4U 1626-67 was 
observed to spin down at a rate v ~ -7 x 10~ 13 Hzs 1 and spin 
up at v w +8.5 x 10~ 13 Hz s 1 . These spin variations correspond 
to torques 

N = 2nvl = 6.3 x 10 25 V-i 3 / 38 Nm, (66) 

where / is the moment of inertia of the neutron star, which we 
measure in 10 38 kgm 2 , and we measure the spin acceleration in 
units of l(T 13 Hzs -1 . The required torque is significantly larger 
than what is produced by N s h em - in any of our models, but iVdyn 
is even an order of magnitude greater than this value if |e| = 0.1. 
This overestimate can be explained by that we have assumed a 
unidirectional fi^dyn across the disc surface, while it might be 
more realistic to expect that the toroidal field is organised in 
magnetic annuli with opposite polarities, or that |e| is signifi- 
cantly smaller fraction of the total turbulent magnetic field. 

We explain the torque reversals in the same way as in 
Torkelsson (1998), where the dynamo undergoes a field rever- 
sal, which in our model corresponds to e changing sign. One 
might speculate that, during this field reversal, the disc passes 
through a state corresponding to e = 0, in which the disc will 
have a smaller inner radius than during the states with an active 
dynamo. The closer to the surface of the neutron star that the ac- 
cretion disc extends, the more of the X-ray emission from the ac- 
cretion columns it can absorb and re-process. Torkelsson (1998) 
showed that the observed time scales on which the torques re- 
main constant in the sources Cen X-3, OAO 1657-415, and GX 
1 +4 are comparable to the global viscous time scales of their ac- 
cretion discs, which constrains the mechanism that is responsible 
for the reversals of the magnetic field. 



viscous stress between the accretion disc and the boundary layer, 
which can transfer angular momentum outwards at a rate that is 
comparable to the one at which it is advected inwards by the 
accreting matter itself. 

We have also found that the dynamo leads to that the inner 
edge of the accretion disc occurs at a radius that is larger than 
the traditional Alfven radius. This effect is weak, though, for a 
realistic value of the dynamo-generated magnetic field. 
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5. Conclusions 

We have investigated the interaction between a magnetic neu- 
tron star and its surrounding accretion disc in the case where the 
accretion disc is supporting an internal dynamo. The magnetic 
field that is produced by the dynamo can lead to a significant 
enhancement of the magnetic torque between the neutron star 
and the accretion disc, compared to what is seen in the model 
by Ghosh & Lamb ( 119791 ). This extra magnetic torque can ex- 
plain the large variations in spin frequency Cen X-3 and OAO 
1657-415 (Bildsten et al. 1997). Furthermore, a reversal of the 
magnetic field that is generated by the dynamo, similar to the 
reversals that we see of the magnetic fields on the Sun, could 
explain the torque reversals in these objects. 

From the way that we calculate the structure of the accretion 
disc, we find two kinds of solutions with different behaviours at 
the inner edge. A few of our solutions have case D boundaries 
at which the density and temperature go to at finite radius, 
while most of our solutions have case V boundaries at which 
the accretion is driven entirely by the magnetic tension between 
the accreting matter and the neutron star. In this case there is a 
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Table 2. The inner edge of the accretion disc and its torque on the fiducial neutron star. 
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